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SUMMARY 


General equations are developed for the design of efficient struc- 
tures protected from thermal environments typical of those encountered 
in boost- glide or atmospheric- reentry conditions. The method is applied 
to insulated heat- sink stressed- skin structures and to internally cooled 
insulated structures. Plates loaded in compression are treated in 
detail. Under limited conditions of plate buckling, high loading, and 
short flight periods, and for aluminum structures only, the weights of 
both configurations are nearly equal. Load parameters are found and are 
similar to those derived in previous investigations for the restricted 
case of a constant equilibrium temperature at the outside surface of the 
insulation. 


INTRODUCTION 


An aircraft which flies at high speeds in, or reenters into, an 
atmosphere generates severe thermal environments by its high-speed 
motion. Elevated temperatures have a deleterious effect upon such 
material properties as the allowable yield stress and Young's modulus; 
and, since the structural strength depends upon these properties, the 
amount of load that the structure can support decreases as the temper- 
ature rises. (in rare instances mild heating might increase an allow- 
able yield stress, but this behavior is not utilized in design.) The 
high temperatures associated with the thermal environment can raise the 
temperature of the aircraft structure to a point where the material has 
insufficient strength to support the aerodynamic loads. Structure tem- 
peratures can be reduced by protecting the structure from its thermal 
environment through the use of insulation or cooling, or both, thus per- 
mitting lighter primary structures. However, the heat- protection system 
itself adds weight to the aircraft and must be considered when the 
weight of the protected structure is calculated. 
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For the purposes of this report, an optimum structure is defined 
as the lightest possible structure which supports a given load within 
a specific thermal environment. In obtaining the structural weight the 
weight of the thermal protection system is included. Hence, a protected 
structure is efficient if the combined weighi of the protected structure 
and protection equipment is less than the weight of an unprotected struc- 
ture designed for the same flight conditions. 

References 1 and 2 consider insulated heat- sink structures wherein 
the structure itself is used as a heat sink, and insulation is placed 
between the high- temperature atmosphere and the structure to retard the 
flow of heat. Throughout those analyses the effective environmental 
temperature was considered constant. It was found that the system 
weight could be expressed as a function of structure temperature. An 
optimum structure was determined by minimizing the weight with respect 
to the temperature. From the results load parameters were derived 
through which it was possible to plot the weights of optimum configure 
tions. The load parameters include the applied load, the flight time, 
and the insulation characteristics. There was a load parameter associ- 
ated with each of the design criteria of yield, buckling, and post- 
buckling failure. The minimization of the weight with respect to the 
temperature works well for the case of constant environmental tempera- 
ture, but involves an excessive amount of algebraic manipulation when 
applied to problems where the thermal environment varies with time. 

A more general thermal environment, which is typical of that expe- 
rienced by a boost- glide missile or atmospheric- reentry body, is con- 
sidered herein. The work of the previous reports is extended to include 
insulated and internally cooled structures as well as the insulated heat- 
sink configurations. The algebraic difficulties of the parametric method 
are circumvented by applying Lagrange ' s method of undetermined multi- 
pliers. The results lead to load parameters very similar to those of 
the case of a constant- temperature environment. Computations are made 
for a water-cooled structure, and the results are compared with those 
for an insulated heat- sink structure on the basis of plate buckling 
strength . 


SYMBOLS 


A constant, °F/(sec) n 

a dummy index 

B exponential decay constant, l/ sec 


b 


plate width, ft 


t-q ON On H 
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Ci heat capacity of insulation, Btu/lb- °F 

C2 heat capacity of metal primary structure, Btu/lb- °F 

Cj heat of vaporization of coolant, Btu/lb 

c* effective heat of vaporization when weight of pump and piping 

^ equipment is proportional to water required, Btu/lb 

E Young’s modulus, lb/sq ft 

e exponential base, 2 ,718 



^(^1^2) func ^i° n (t-jJ and (^ 2 ) 


H = BG 


h 

K 

k 

n 

P 

Q 

T 

t 

u,x,y 

W 

w 


(see eq. (13)) 

heat- transfer coefficient, Btu/sq ft-sec-°F 

buckling coefficient 

thermal conductivity, Btu/ft-sec-°F 

constant 

load, lb 

amount of heat absorbed by coolant, Btu/sq ft 
function of (t^ and (t*^ 

temperature, °F 
thickness, ft 

dummy variables of integration 
weight per unit area, lb/sq ft 

total weight per unit area of pump and piping equipment when 
distributed over heated surface of vehicle, lb/sq ft 
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V weight per unit area of fixed weight portion of pump and 

P piping equipment when distributed over heated surface of 

vehicle, lb/sq ft 


a 


c 2P2 t 2 t l’ 


l/ sec 


P 

7 


€ 


A 

P 

P 

a 


T 


cr 


* 


% 


e 


S’, 


e 


exponential decay constant (associated with structural mate- 
rial properties), l/°F 

constant (associated with structural material properties), 
lb/sq ft 

ratio of variable portion of pump and piping weight to cool- 
ant weight (see eq. ( 52 )) 

nondlmensional time parameter, <xt 

Poisson's ratio 

density, lb/cu ft 

general material strength property or parameter, lb/sq ft 
yield stress, lb/sq ft 


time, sec 

critical time, sec 


Lagrange undetermined multiplier 


load parameter for buckling, 

lb \ 2//? / Btu V 1 /lb V 1 
ft 2 / \ft3-°F/ \ft 5 / 



Bb 


2 


k lPl' 


load parameter for yield. 



f Btu \ 

,ft3-°Fi 




lb 


2 


k lPl’ 


lb 

ft 3, 


-l 
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Subscripts: 

1 insulation 

2 primary structure 

3 coolant 

aw adiabatic wall (temperature) 

b buckling 

eq equilibrium temperature 

i indie ial notation 

j.,k signify beginning and end of coolant vaporization period 

max maximum 

o initial conditions 

v vaporization 

y yield 


ASSUMPTIONS 


A sketch of the insulated heat- sink stressed- skin structure is 
shown in figure 1. Figure 2 is a sketch of the insulated and cooled 
structure. The primary structures are considered to be edge- supported 
plates loaded in compression such as covers of a box beam, or some other 
load- carrying surface exposed to aerodynamic heating. The volume of 
coolant required per unit of heated surface area is denoted by a 
"thickness" ty The subscripts 1, 2, and 3 have been assigned to 

the insulation, primary structure, and coolant, respectively. 

The following assumptions have been made to simplify the analysis: 

(1) The heating is uniform over the entire surface of the member j 
that is, the heat flow is one dimensional. 

(2) The primary structure supports the entire load. 

(3) The temperature gradient through the metal skin is negligibly 
small . 
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(4) For the case of the heat- sink structures the thermal capacity 
of the insulation can be neglected with respect to the thermal capacity 
of the primary structure (metal skin). This approximation may be 
improved by the iteration method of appendix A of reference 1. The 
iteration method is valid for cases where CjjD^t^ < 2c2P2^2 should 


be used whenever 


c ipi t i > _i_ 
c 2 p 2 t 2 20 


For the internally cooled structure it 


is assumed that all of the heat passing through the outside surface of 
the insulation is absorbed by the coolant. 


(5) The outside surface temperature of the insulation T e( ^ is a 

known function of time. This surface temperature can be determined by 
the method of appendix B of reference 1. 


(6) The heat capacity and thermal conductivity of the materials 
are independent of temperature. (Note that the yield stress and Young’s 
modulus are considered functions of temperature.) 

( 7 ) The coolant used in the insulated and cooled structure absorbs 
heat by vaporizing. The thermal capacity obtained by raising the cool- 
ant temperature from 75° F to 212° F is neglected. 

Assumption (4) has the greatest effect upon the accuracy of the 
analysis of the insulated heat- sink configuration. For the ranges of 
thermal capacities commonly encountered in practice, the assumption will 
usually lead to acceptable engineering accuracy without iteration. This 
accuracy is evaluated in reference 3- 

Assumption ( 7 ) will lead to a conservative (high) answer for the 
weight of the insulated and cooled configuration. The tendency to over- 
estimate the coolant required is somewhat offset by considering the 
heating period to exist only during the period wherein the equilibrium 
temperature exceeds the coolant boiling temperature. This situation is 
discussed subsequently in more detail. 


CHARACTERISTICS OF THE ASSUMED THERMAL ENVIRONMENT 


References 1 and 2 examine the problem :>f finding the optimum 
weight of an insulated structure exposed to a constant equilibrium tem- 
perature at the outside surface of the insulation. In order to optimize 
the weight of a structure exposed to boost- glide or reentry flight paths, 
it is necessary to consider a variable equilibrium temperature at the 
outer surface of the insulation. The development of the equations in 
parts of this report is general and independent of any particular time 


h vo vo ir 1 
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variation of the equilibrium temperature except for the unobstructive 
requirement that the equilibrium temperature eventually decreases during 
the flight period, and that, consequently, the primary structure temper- 
ature has a definite maximum value. This mathematical restriction does 
not impede the solution of physically real problems concerning recover- 
able vehicles which eventually must slow to landing speeds and hence, as 
a natural result, experience a decrease in surface equilibrium temperature. 

In order that demonstrative results might be obtained to illustrate 
the application of the optimizing relationships, the following particular 
time variation in equilibrium temperature T e ^ was chosen: 

T eq ' T o = AT n e~ BT (l) 

The constants A, B, and n may be determined to give a best fit to 
some experimental or predicted surface- temperature curve of interest. 
Relations which may be useful for curve fitting are 

Bt = n (2) 

when 


and 



or 



(5b) 


Equation (l) is plotted in figure 3 for five values of n to show how 
the shape of the curve is modified by n. In addition, figure 4 shows 
by a nondimensional plot that the curve has a sharper peak as n is 
increased. 


This particular temperature variation was chosen for a number of 
reasons: As seen in figure 3, the equation generates a temperature 

history which might be considered typical of a boost-glide or reentry 
heating condition. In addition, for the purposes of analysis it was 
desired to have an analytic function which was easily differentiable 
and integrable . Finally, the number of constants required to fit the 
curve was kept small thereby keeping the constants appearing in the 
loading and weight parameters small. 
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Reference 3 gives several solutions to the heat- transfer problem 
of the insulated heat- sink structure. For the purposes of this investi- 
gation the first approximation, as was mentioned previously in the 
assumptions, was used. Equation (3^) of reference 3 gives the primary 
structure temperature as 


where 


and 



A = 


k l T 


c 2 ^ 2 ^ 2^1 


5 - 


k i 

ht*L 


It is shown in reference 3 that this solution for aerodynamic heating, 
based on adiabatic wall temperature and heat- transfer coefficient, can 
be converted to the solution of the problem of known surface tempera- 
ture by performing the manipulation that, as h - >00 > T aw^ T eq‘ See 

ref. 1 or 3 for the determination of T eq . This manipulation results 
in 

t 2 “ T o = ( T eq - T o )^_ 0 (l - e“*) + [l - e - ( A-u)j d ( ^~ T o) du (4) 


Equation (k) may be written as 


T 2 - ' (V - + / q T t 


dy 


( 5 ) 


where 


a = 


c 2 P 2 t 2 t l 


Through the definition of a equation ( 5 ) shws that the first approxi- 
mation will always lead to the general result that, functional! y 


H vo MD It 1 
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T 2 = T 2 (t 1 t 2 ) 


( 6 ) 


where 


(tit 2 ) always appears as a product. 


Because of this characteristic, 


2 Uai ' HatJ 


where the subscripts on the parentheses indicate quantities held con- 
stant during differentiation. Use will be made of this property 
subsequently. 

Combining equation (l) with equation ( 5 ) yields 


t 2 - 



( 8 ) 


Temperatures obtained from equation (8), with n = 2 and 

^eq,max = 3,000° F, are shown as a dashed line in figure 5- For the 

case of the heat- sink configuration the maximum temperature of the 
primary structure occurs at the time when its temperature equals the 
equilibrium temperature; thereafter, aerodynamic cooling takes place. 
The maximum temperature of the primary structure is the design temper- 
ature and may be found by setting the right-hand side of equation (8) 
equal to the right-hand side of equation (l). The result is 

x( T 2,max - T 0 ) - (BT cr ) n =‘ ETcr (9) 

where BT cr is determined from 

l(Br cr ) n = / Tcr y^e^Ky (10) 

In equations ( 9 ) and (10) the time of occurrence of the maximum 

temperature T P is denoted as the critical time t . Mathemati- 
cs, max cr 

cally there are three unknowns in the two equations, T cr’ 

and a. The value of a depends upon the insulation and primary- 
structure thicknesses ti and t2 which are not known as yet. The 

additional equation required to determine the proper value of a is 
found from the condition that the structure must be an optimum. Equa- 
tions (9), (l0), and the optimization equation are sufficient to deter- 
mine the unknowns. 
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In the following development the weight of an optimum heat- sink 
structure is compared with that of an optimum cooled structure. It is 
therefore desirable to be able to express the weight of the cooled con- 
figuration in terms of the same parameters as those that are obtained 
for the heat- sink configuration, and the following equations will be 
useful for the manipulations of the equation for the weight of the 
cooled configuration. 

A result of assumption (7) is that the coolant is considered to be 
expended only during the period in which the equilibrium temperature 
exceeds the temperature at which the coolant vaporizes. The times at 
which the equilibrium temperature equals the vaporization temperature 
may be found from 

TV - To - T eq ' T o * ^< Bl (ll) 

The vaporization temperature and initial temperature are known as design 
conditions, and the two roots of equation (ID.) are found as Btj and 

BT k . (See fig. 5.) These roots are functions only of A/B n , a known 

characteristic of the equilibrium- temperature function used herein. 

The amount of heat absorbed by the coolant in 



Because t 1 is not known until the optimum amount of insulation has 
been determined, it is convenient to define G as follows: 



in the general case where G is proportiona. to the amount of required 
coolant. For the particular equilibrium- temperature variation con- 
sidered herein, because of equation (ll), the general equation (12a) is 
better expressed as 



HVOMD f 
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or 


G 


H 

B 


where 


H 






dx 




Br k + 


E 

0=1 


n! 

1 

(n - a) ! 

k BT j) a 



( 13 ) 


The second expression is obtained by successive integration by parts of 
the first expression. 

The values of Btj and Bt^ are functions only of A/B n as was 

noted after equation (ll), and thus H is a function only of A/B n . 
Because of equation (3b), A/B 11 is a function only of T e ^ max - T Q 

and n, and it is assumed that T v , 1*eq,max> T 0 , n are 

known. 


INSULATED HEAT- SINK STRUCTURES 


General Equation 

The general case of the insulated heat- sink stressed- skin structure 
subjected to variable temperature at the outside insulation surface (see 
fig. l) is discussed first. The insulated panel is heated aerodynami- 
cally, and part of the heat entering the insulation outer surface from 
the boundary layer is radiated away to the surroundings. The remainder 
of the heat input is conducted through the insulation to the structure. 
Appendix B of reference 1 discusses a simple method of determining the 
temperature of the insulation outer surface by considering a heat balance 
between input and radiation. This surface temperature is called the 
equilibrium temperature. The problem is to find the combination of 
insulation and stressed- skin thicknesses which results in the lowest 
combined weight for a given loading and temperature environment. Because 
the weights of surfaces are concerned, it is only necessary to determine 
the minimum weight per square foot of surface to find the optimum dimen- 
sions of the insulated panel. The equation relating the weight to the 
thicknesses of the insulation and the primary structure is 
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W - + P2^2 

which may be written functionally as 

g(ti,t 2 ) = 0 = W - p 1 t 1 - o 2 t 2 (l^) 

The assumption is made at this point that the kind of metal to be 
used in the primary structure and the kind of insulating material have 
been chosen. (The insulating material should oe selected on the basis 
of minimum k-^p^ and ability to withstand T e ^ max as is discussed in 

refs. 1 and 2.) The following physical and thermal characteristics are 
thus fixed: k^, pj_, P 2 , and eg* It is also assumed that the dimen- 

sion b is prescribed. (See fig. 1.) 

The allowable load depends upon one or more material properties 
and the thickness of the primary structure; therefore, 

P = P(c,t 2 ) (15) 

where P is the imposed load for which the plate must be designed. The 
load and plate width b are prescribed in a design problem, and T 2 is 

usually adjusted to make the design load equal the allowable load. The 
symbol a stands for a general material strength property such as 
Young's modulus or yield stress, and is assumed to be a function of 
tempe r atur e alone. Thus , 

a = o(T 2 ) (16) 

The temperature T 2 is that of the primary structure. For design 
purposes the value of a is determined at the maximum value of T 2 

because the metal mechanical properties are at a minimum at the maximum 
temperature of the primary structure. The maximum value of T 2 is not 

yet known because the relative dimensions of the optimum structure are 
yet to be determined. The temperature of the primary structure T 2 is 

a function of time and the insulation and metal thicknesses; thus, 

T 2 = T 2 (t x , t 2 , t) (17) 

The temperature T 2 is also a function of the time- dependent equilibrium 
temperature. However, it is assumed that T e ^ = is a given func- 

tion, and therefore T e(1 does not appear explicitly in equation (17). 

In references 1 and 2 the determination cf the optimum weight was 
approached as an extremum problem where the weight was to be minimized 


CT\ CT\ H 
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for a given imposed, load. Identical results may be obtained by maxi- 
mizing the load with respect to a fixed weight. Lagrange's method of 
undetermined multipliers is used herein, with the weight considered 
fixed and the load to be maximized. The thicknesses t^ and tQ are 

considered to be the variables, and equation (lU) is used as the equa- 
tion of constraint. The following two equations result from the appli- 
cation of Lagrange's method to equations (l4) and (15): 




0 


(l8a) 


+ * = 0 
0^2 0^2 


(18b) 


where \|r is the Lagrangian multiplier. The set of equations (l8) is 
subject to the restraining condition of equation (l4). Equations (l8) 
become 


which yield 



= 0 



*p 2 = 0 


P2 

Pi 


dP/dt 2 

dP/dt-L 


( 19 ) 


The partial derivatives in equation (19) can be expanded by using 
the previous functional relationships with the objective of obtaining 
partial differential terms that may be easily evaluated from known 
strength and heat- transfer equations. From equation (15) (keeping in 
mind eqs. (l6) and (17) and considering t^, t2, and r as the inde- 

pendent variables). 




r 


where the subscripts indicate variables held constant during the 
differentiation . 



Ik 


It is true that 

2-\ 


So 


da 




\ 3t 2/ tl ,T “aV^a/t^T iT 2^2 / tl( , 

where the last form is a consequence of equation (l6). Thus, 


dP 

dt 2 


fbP 

i dtc 


bp\ da [^2 


[bo 


dTp\< 


(S) 


t X ,T 


which is a convenient form. 
Similarly, 

dP . 
Sti 


SP 


da 




\bo/+ dTo\St n 


2 J -'t 2 ,T 

Substitution of equations (20) into equation ( .9) yields 

^dP\ da 2 s 


P 2 


< bp > 
idt 2 y 


Pl ($E\ dCT /dTg\ 

H dT2 ^i 4 2 ,t 


^a/t dT 2 \dt 2 / 

C Uq, T 

ctf>\ da /dTg\ 


Wt 2 (iT 2\ at l/ t ^ T 


which may be reduced to 


r bP > 
v St 2y 


'&V 


ti, T 


Pl /dP\ da /^T 2 


sr/ t2 dT 2 Vd tl / t2jT \S a / to . T 


&|V 


u 2 , 


Equation (21a) is multiplied by t 2 ^t]_ to obtain 
^2^2 


?/ap \ 
Uto / 

t te) 

2 ^v 

da fb T 2 \ 

4%,, 

+ / dT?\ 

‘ - L/t 2> T 


(20a) 


(20b) 


(21a) 


(21b) 
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Use will now be made of the property of the heat- transfer equation 
exhibited by equation (7)- When equation (j) is substituted into equa- 
tion (21b) there results 


P2^2 

01*1 


*2 /dP \ 

ti \dt 2 j 


dP\ do /dT 2\ 

dcr/ + dtoldt-, / 

/t 2 2 \ l/t 2 , T 


+ 1 


( 22 ) 


It is only necessary to substitute the appropriate partial deriva- 
tives into equation (22) to obtain the optimum conditions for an insu- 
lated heat-sink structure. Equation (22) is a general result and applies 
to any time variation of equilibrium temperature. Examination of the 

first term on the right-hand side of the equation shows that < 0 


because material properties generally decrease with temperature, and 
/ dT 2 ' t 


that 


< 0 because an increase in the insulation thickness will 


decrease the temperature T 2 - All other terms are positive. Equa- 
tion (22) thus indicates that 


p 2^2 ^ ^ 1^1 


or that the primary structure is heavier than the insulation for an 
optimum heat- sink design. 


Application of Mechanical Strength Criteria 

Design on the basis of yield stress .- The optimization equation 
(eq. (22) J is used to detenSne the opt imum configuration for a struc- 
ture designed on the basis of yield. For the yield condition equa- 
tion (l5) becomes 


P = bt 2 a y 


(23) 


An empirical relationship that relates the yield stress to temperature 
was used in references 1 and 2 and will be used herein^ the specific 
form of equation (l6) to be used in the analysis is then 



(2k) 
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This empirical representation of the yield stress is plotted in figure 6 
with experimental data for 2024- T3 aluminum alloy. The functional equa- 
tion ( 17 ) Tor the structure temperature is represented by the specific 
equation (8). For design purposes, however, the maximum structure tem- 
perature should be used. (This maximum is given by the simultaneous 
solution of eqs. ( 9 ) and (10).) Taking the appropriate partial deriva- 
tives of equations ( 9 ), (10), ( 23 ), and (24), and substituting the 
results into equation ( 19 ) yields 


B 

a 





( 25 ) 


where the subscript e is used to distinguish the load parameter 


a 





from the load parameter 



b 


2 


klPlT* 


which is 


a result of references 1 and 2 where the equilibrium temperature was 
considered constant. The two parameters are the same dimensionally. 

In the case of constant equilibrium temperature the flight period (and 
hence the heating period) is known. In the variable equilibrium temper- 
ature considered herein the time of occurrence' of the maximum structure 
temperature is not known a priori, but the value of B is known from 
the characteristics of the equilibrium- temperature curve. 


The value of a for the optimum configuration can be found from 
the simulataneous solution of equations ( 9 ); <10), (24), and ( 25 ) (with 
T 2 max used in eq. (24)). For the purpose of computing design charts 

it is easier to solve this system of equations inversely by first 

assuming values of B/a for fixed values of a/b 11 , determining the 
values of T 2 max and cfy from equations ( 9 ), (10), and (24) (or plots 

thereof), and then computing Trom equation ( 25 ). After these 

final results are plotted a may be determined for given values of Qy e * 


The substitution of equations ( 23 ) and ( 25 ) into equation (l4) yields, 
after manipulation 


/Bb 2 ^ 

1/2 

-£2 

l + 

( a y) 

2 

B 1 

\ k lPl) 

1 °y 

\P2j 

a c 2 n y, e 


1/2 


y,e 


(26) 


where Oy and B/a are found through charts of the solution of equa- 
tion ( 25 ) by using the given value of These values of a v and 

y> e y 

B/a are values for an optimum structure because of equation ( 25 ). 
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Design on the basis of plate buckling .- For plate buckling the equa 
tion relating the applied load to the plate dimensions and material pro- 
perties characterized by equation ( 15 ) is 


P 


Krt 2 E V 

12 (l - H 2 ) b 


(27) 


where K = 4 for a plate simply supported along the edges. (See ref. 4 
An empirical representation for the elastic modulus is (from ref. 2) 


E 


= 7 


i,b c 


-Pi,b T 2 


( 28 ) 


where 7± f \ > Pf b are a lj us t e( l to give good fit for two ranges of 

temperature. The representation is shown in figure 7 for aluminum. 


Equations ( 9 ), (10), ( 27 ), and ( 28 ) are substituted into equa- 
tion (22) to obtain the relationship between the buckling load factor 
and the optimum value of a. Thus, 



= <V (29) 


It should be noted that e f- e because the loading index P/b^ 
is raised to a different power. 


Substituting equation ( 29 ) into equation (l4) yields, after 
manipulation, 


W/ Bb 2 \ 

b \ k lPl/ 


1/2 




r 

« 2 e 

2/3) 

p 2 


1 . B . 

^(l - 

( 

jt 2 E 

_3 (l - n 2 )_ 

1/3 

a 

V 

P2 2c 2 fi b,e 

'J 




1/2 


(30) 


which is similar in form to equation (26). 
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INSULATED AND COOLED STRUCTURE 


General Equation 

For the development of an optimization equation for the insulated 
and cooled structure, the configuration is assumed to be as represented 
in figure 2. The ” thickness" t^ represents the volume of coolant per 

square foot of heated structure surface and is used with the coolant 
density to express the weight of the coolant. By determining the thick- 
ness of the coolant layer for an optimum configuration the proper weight 
of the coolant can be found. The weight per square foot of insulation, 
structure, and coolant is 

W = p 1 t^ + p 2 t 2 + p^tj + V (31) 


where w is the distributed weight per square foot of outside surface 
of the pump and piping equipment necessary tc circulate the coolant. 

The pump and piping weight is assumed to vary linearly with the coolant 
weight; thus, 


V = Vp + (32) 

where Wp is a fixed weight and e is a constant of proportionality. 

The substitution of equation (32) into equation (31) results in 

W = P^t-^ + P2^2 + ^3^3^ + € ) + w p (33) 

The coolant is assumed to protect the lead-carrying structure by 
absorbing the heat transmitted through the structure by vaporization of 
the liquid. The result of this assumption is that the maximum structure 
temperature (which is taken as the design temperature) is the boiling 
temperature of the coolant. In this case, the problem is to minimize 
the weight for a fixed load, with respect to the variables, t^, t 2 , 

and tj. 

If the thermal capacity of the insulation and the structure is 
small compared to the latent heat of vaporization of the coolant, the 
heat- transfer problem is essentially one of quasi- steady state with all 
the heat passing through the insulation being absorbed by the coolant. 
When T e( ^ < T y no coolant is boiled. The conservative assumption is 

made that heat is absorbed only by the vaporization of the coolant. 

Thus, the heat- transfer equation may be written as 


t-q C\ ON H 
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& / k ^ " Tv ) dT = c 3 p 3 t 3 


(3*0 


T j 

where Tj = t % is the period during which T e( ^ = T y . 
Equation (3M may be represented functionally as 

/ T k 

( T eq ” T v) dr 

r J 


(35) 


Equation (35) is the equation of constraint which is used in conjunction 
with equation (33) to minimize the weight. 

When the Lagrange multiplier technique is applied to equations (33) 
and (35) there results 

dW dq 

-r + + = 0 

dtp dtp 


3L + * sa. . 0 

^t2 St2 


|«- + *| 2-.0 

dtj dtj 

From equation (35) it is found that 

dq 


dtr 


= 0 


and therefore 


dW_ 

dtc 


= 0 


or that the extremum problem is independent of t 2 - This result is to 

be expected because t 2 is in actuality dependent only upon the struc- 
ture temperature, the imposed load, and the value of the temperature- 
dependent material property in the strength equation, all of which are 
prescribed and therefore invariant. 
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The remaining two equations yield 


Sw/St]^ dq/dti 

cW/dtj dq/dtj 


( 56 ) 


Taking the appropriate derivatives in equations (33) and. (35) and 
substituting into equation (36) yields 


P^_t]_ ~ P^tj(l + e) 


( 37 ) 


which is the optimizing relationship. 

The weight per square foot of surface area becomes 


W = 2p 1 t 1 + p 2 t 2 + w p 
— p 2 t 2 4* 2 


^ipiC 1 + 0® 


+ w_ 


(38) 


where 


G = 



- T v )dT from equation (l2a). Equation (58) 


represents the weight of an optimum structure subjected to any arbitrary 
variation in equilibrium temperature and is not restricted to the partic- 
ular equilibrium- temperature variation used herein. Equation (58) also 
shows that the proportionality constant between the pump and piping- 
equipment weight and the coolant weight effectively reduces the thermal 
capacity of the coolant to a value which is 


C 3 

3 ” ITT 


(39) 


Application of Mechanical Strength Criteria 

Design on the basis of yield stress .- The optimum weight of a 
structure designed on the basis of either tensile or compressive yield 
will be determined. Equation (38) shows thet the weight is linearly 
dependent upon the thickness of the primary structure. From equation (23) 


t 2 
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and therefore 


W 


= 


= P2£ + 2 J k iPi(i + £ ) g 


bO, ; 


(40) 


In equation (40) a is evaluated at the boiling point or evaporation 

v 

temperature of the coolant. It should be noted that equation (40) is 
not restricted to a particular mode of time variation of T e( ^ if G 

is not so restricted. 


The weight of an insulated heat- sink structure may be compared to 
that of an insulated and cooled structure. Equation (40) can be 
expressed in terms of the load factor Q as 



V J 


Since H is a function only of T e( ^ max and n, the weight parameter 
can be plotted against T e q,max ^y e as the case 0 f the heat- 

sink configuration. Since the left-hand sides of equations (26) and (4l) 
differ slightly because of Wp, minor arithmetical computation is neces- 
sary to compare weights directly unless Wp = 0. If Wp = 0, equa- 
tions ( 26 ) and (4l) can be compared directly. 

Design on the basis of plate buckling .- The application of equa- 
tion ( 27 ) (with K = 4) to the buckling criterion follows in a manner 
similar to that which was used for the heat-sink structure. From 
equation ( 27 ) 



Substitution of equation (42) into equation ( 38 ) yields 


W 


w p = bp 2 


/ p\ 3(1 - n 2 ) 

1/3 f 

X O 1 

\b 2 / rt 2 E 

1 


kqPl (l + e)G 


(42) 


(^5) 


Equation (43) is manipulated to compare with the insulated heat- sink 
structure expressed by equation (30)* Equation (43) then becomes 
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and, as before, the weight per square foot of structure can be com- 
pared directly by comparing values of the weLght parameters when 



COMPUTATIONS 


Structural weight and loading parameters have been computed which 
permit a limited comparison to be made between an insulated heat- sink 
structure and an insulated and cooled structure. The computations were 
made to determine whether or not significant differences exist between 
the efficiencies of the two configurations. 

In order to obtain comparative results it was necessary to assign 
certain characteristics and values to the parameters which determine 
the configuration weights. Equilibrium- temperature variations corre- 
sponding to boost- glide of reentry flight paths as given in equation (l) 
were used with n = 2. This temperature variation with time is shown 

in figure 5- Various values of the ratio a/b^ were chosen so as to 
vary the maximum equilibrium temperature. (See eqs. (5)-) 

The parameters for the insulated and cooled structure were con- 
sidered first. It was necessary to choose a cooling fluid. Several 
coolants were considered briefly, and water was chosen on the basis of 
availability, simplicity of handling equipment, and high latent heat of 
vaporization. The water was assumed to boil at atmospheric pressure to 
limit the design pressure of the coolant passage to one atmosphere. 

Thus, the coolant boiling temperature was fixed at 212° F. For this 
operating temperature 2024- T3 aluminum alloy was chosen as a familiar 
and practical structural material. Somewhat arbitrarily, the propor- 
tionality constant between the pump and pipiig weight and the coolant 
weight was taken as 


e = 0.072 

(which is the same as taking c^ = 900 Btu/lo in eq. (39))- 


H vo V Of 1 
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With these chosen values equation (44) was used to compute the 
weight parameter for the insulated and cooled structure for the buckling 
criterion. The results are shown in figure 8 where the weight parameter 
is plotted as a function of T e( ^ max and e* 

The 2024- T3 aluminum alloy was also chosen for the structure mate- 
rial in the heat- sink case even though, if efficiency were the only con- 
sideration, the results of reference 2 indicate that HK31A magnesium 
alloy would be better. Aluminum was chosen to permit a more direct com- 
parison between the heat- sink and cooled configurations (and because, in 
considering present fabrication practices, aluminum seems to be favored 
over magnesium for structural members near the external surface of a 
heated aircraft). Again equilibrium- temperature variations given by 
equation (l) with n = 2 were used. For the heat- sink structure suc- 
cessive values of the maximum equilibrium temperature (and hence A/B^) 
were chosen and then, for each value, a range of B/a was assumed. 

The critical time and the maximum structure temperature were determined 
from equations (9) and (10). The load parameter was then computed from 
equation (29) and is plotted as a function of maximum equilibrium tem- 
perature and B/a in figure 9- This figure, together with figure 7 
and equations (9) and (10) to determine the maximum structure tempera- 
ture, was used with equation (30) to compute the weight of the insulated 
heat-sink configuration. This weight is plotted in figure 10 in the 
form of the weight parameter as a function of T e( ^ max and 

In addition to the previous computations the effect of varying the 
boiling point of the water was examined. Ttie necessary increase in the 
weight of the coolant channel with an increase in pressure was neglected 
In general, it is better to boil water at the highest possible tempera- 
ture to obtain a slightly larger total heat absorption and to reduce the 
temperature gradient through the insulation. This gain was slight, how- 
ever, and probably would not offset the increased weight of the coolant 
channels in an actual design. It should be mentioned that when other 
factors such as the internal environment of the space vehicle are con- 
sidered the best overall efficiency may be obtained by vaporizing the 
water at room temperature (75° F). Although this may not be optimum 
from the strict structural- efficiency viewpoint, operation at this low 
temperature eliminates the need for internal cooling equipment to pro- 
vide a habitable temperature. 


RESULTS AND DISCUSSION 


General 

The relative weights of the insulated heat sink and the insulated 
and cooled configurations with aluminum used as the structural material 
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may be observed by comparing figure 10 with f:.gure 8, This comparison 
is limited to cooled structures where the wei<-£ht of the pumps and piping 
may be considered proportional to the coolant weight, with a factor of 
proportionality of 0.072. For these calculations there was very little 
difference between the weights of the two configurations; in fact, for 
the range of e considered for the heat- sink structure, figure 10 

may be superposed upon figure 8. As decreases the insulated and 

cooled configuration becomes more efficient. 

It is worth noting that for the insulated heat- sink configuration 
the weight of the primary structure is always greater than the weight 
of the insulation in an optimum design. This fact is shown by the par- 
tial differential equation (eq. 22) which defines the necessary relation- 
ships for an optimum structure. The condition which determines the 
optimum for the insulated and cooled configuration is that the weight 
of the insulation is equal to the weight of the coolant when the coolant 
weight is calculated by using the reduced the:mial capacity. (See 
eqs. ( 52 ), (57), and (59)-) 

There is a significant difference between the relationships that 
exist for the optimum heat- sink configuration and those for the optimum 
internally cooled configuration. The equations which determine the 
optimum insulated heat-sink structure (eqs. ( 25 ) and ( 29 )) contain the 
applied load P, whereas the equation for the internally cooled struc- 
ture (eq. (37)) is independent of the load. In the former case the 
amount of heat sink depends upon the plate thickness, and, because the 
heavier loads require thicker plates, an increased amount of heat sink 
is more readily available. For the internally cooled configuration the 
structure temperature was prescribed so that ihe amount of absorbable 
heat became independent of the thickness of the primary structure. Thus, 
the optimization equation (eq. (37)) relates only the insulation weight 
to the coolant weight and is independent of the applied load. The sim- 
plification introduced by prescribing the structure temperature is 
impossible for configurations such as the heai-sink configuration if a 
true optimum is sought. 

A similar but more subtle point arises for internally cooled con- 
figurations when several different coolant vaporization temperatures 
are under consideration. In general, the coolants will have different, 
but prescribable, vaporization temperatures. Because the structure 
thickness t 2 depends upon the temperature- dependent material prop- 
erties, t 2 will vary with the coolant vapor:. zation temperature chosen. 

The problem of finding the coolant which provides the least total weight 
of structure now depends upon t^. However, :.f each coolant and its 

vaporization temperature is considered independently the equations devel- 
oped herein are applicable, and the optimum weights for each coolant can 
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be found. A comparison of these optimum weights (by directly comparing 
the values obtained for the weight parameters) will reveal the most 
efficient coolant. 

As was mentioned in the section entitled Computations, the slight 
increase in themodynamic efficiency realized by boiling the coolant 
(water) at a pressure higher than atmospheric pressure would probably 
be offset by the increase in the weight of the coolant passage required 
to contain the higher pressures, and, if overall vehicle efficiency be 
considered, the total vehicle weight might be minimized by vaporizing 
the water at 75° F to eliminate the need for heavy internal climate- 
conditioning equipment. 


Use of the Computed Curves 

All of the computed curves contained herein are based upon the 
thermal environment described by equation (l), with n = 2. As was 
mentioned previously, curves for other heating conditions may be com- 
puted by the procedure outlined at the end of the section entitled 
Computations . 

Once a set of curves has been obtained for the desired value of n, 
the optimum protection system may be found for any combination of A 
and B. The first design step is to determine the type of insulation to 
be used. The optimum insulating material will be that which, first of 
all, will withstand the peak value of the equilibrium temperature 
T e q, max j and second, will have the minimum value of k^p^. It is obvious 

that this insulation will be an optimum for both the heat sink and the 
internally cooled configurations, whichever may later turn out to be 
the better method of designing the structure. The values of B, b, P, 
and kppp are used to compute the loading parameter or ft-u. , 

whichever may apply ). The structural configuration of least weight 
(heat sink or internally cooled) may be obtained from a comparison 
of the weights determined from charts such as figures 8 and 10. 

Under the assumption that the weight comparison, together with other 
considerations (such as simplicity of design) indicates a choice of the 
heat- sink configuration, it is then necessary to find the values of t^_ 

and t2 for a minimum- weight structure. The value of B/a for optimum 

conditions can be found from a figure similar to figure 9* Since B is 
known, a is determined; BT cr may be found from equation (10); and the 
maximum temperature of the primary structure Tg max can be determined 

from equation ( 9 ). Once T 2 max is known, either Oy or E (which- 
ever is applicable) is known; and, since the design load is prescribed, 
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the thickness of the primary structure t 2 is determined from either 

equation (25) or equation (27). Because a = — , t-, can easily 

c 2P2^2^1 

be computed. 

On the other hand, if an internally cooled structure is used, 
either Oy or E is Immediately fixed by the coolant vaporization 

point, and t 2 may be readily computed. The times at which the vapor- 
ization of the coolant commences and ends are found from equation (ll). 
By combining equation (5*0 with equation (57)> 


from which 

Cx(l + c) 2 f T k 

k^i ( P 5 t 5) J t ( T eq - T v) dT 

may be computed. This value is used to find t^ through 

equation (37)* 



CONCLUDING REMARKS 


The Lagrange multiplier technique has been used to determine the 
design conditions for optimum insulated heat-sink structures and for 
optimum insulated and cooled structures. This method is advantageous 
over the direct method used in the references when the equilibrium tem- 
perature at the outside surface of the insulation is other than a very 
simple function of time. 

Computations were made to illustrate the application of the analysis 
to a thermodynamic flight path similar to that of atmospheric reentry. 

The computations were restricted to aluminum load-carrying structures. 
Except under conditions favorable to the insulated heat- sink design (high 
structure loads and short heating periods) where the weights of the two 
configurations are equal, water-cooled structures are more efficient theui 
heat- sink structures . 

The results of the calculations are presented in graphic form by 
means of load parameters which account for the loading condition, heating 
condition, and insulation properties. The load parameters are similar 
to those derived in previous investigations for the simpler case of con- 
stant temperature at the outside surface of the insulation. 

The differential equations show the general result that, for optimum 
heat- sink structures, the weight of the primary structure always exceeds 


HVO^O t T 1 



H* MD M3 tT 1 


27 


the insulation weight. For the internally cooled structure an optimum 
configuration exists when the insulation weight equals the combined 
weight of the coolant, pump, and piping if the pump and piping weight 
is considered to vary in linear proportion with the amount of coolant 
required; the weight of the primary structure does not influence this 
optimizing relationship. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Air Force Base, Va., December 13, 1961 . 
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Figure 5*- Plot of equilibrium temperature, pr ima ry- structure temperature, and coolant boiling 

temperature for n = 2 and T n = 3,000° F. 
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Figure J.- Modulus of elasticity as a function of temperature for 
2024- T 3 al umi num . Solid line expresses the analytic approxima 
tion to experimental data. 
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Figure 9-- Load parameter H-j^e for insulated heat- sink aluminum structures plotted against 

■^eq > max ®/ a f° r s-n initial temperature of 75° F. 
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Figure 10.- Weight parameter for optimum insulated aluminum plates designed on the basis of 
buckling plotted against T e q,max an( i ^b e f° r 811 initial temperature of 75° F. 




